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The scale locality of energy fluxes for magnetohydrodynamics (MHD) is investigated numerically for station- 
ary states of turbulence. Two types of forces are used to drive turbulence, a kinetic force that acts only on the 
velocity field and a kinetic-inductive forcing mechanism, which acts on the velocity and magnetic fields alike. 
The analysis is performed in spectral space, which is decomposed into a series of shells following a power law 
for the boundaries. The triadic transfers occurring among these shells are computed and the fluxes and locality 
functions are recovered by partial summation over the relevant shells. Employing Kraichnan locality functions, 
values of 1/3 and 2/3 for the scaling exponents of the four MHD energy fluxes are found. These values are 
smaller compared with the value of 4/3 found for hydrodynamic turbulence. To better understand these results, 
an in depth analysis is performed on the total energy flux. 



PACS numbers: 47.65.-d. 47.27.ek, 47.27.er 

I. INTRODUCTION 

Energy transfer functions in turbulence are the result of 
nonlinear interactions among different scales of motion. Al- 
though all the scales of the flow are coupled together in the 
same manner, through the nonlinear terms, the global trans- 
fer of energy to a scale is dominated by contributions from 
particular scales. Placing adequate bounds on the scales that 
bring the most contributions to the energy flux through a scale 
is crucial for the developing of adequate turbulence models, 
like Large Eddy Simulations (LES), Shell Models and for ad- 
vancing physical understanding of turbulent phenomenology. 
Finding these limits is the main object of locality analysis for 
strongly coupled nonlinear systems. 

In general, in an effort to understand the behavior of scale 
coupling for any system, the resulting transfers due to the non- 
linear terms are investigated for global conserved quantities. 
The redistribution nature of such conserved quantities in spec- 
tral space provides insights into the nonlinear dynamics of the 
system. In the absence of such global conserved quantities, as 
in the case of dissipative systems, ideal invariant quantities are 
used instead, i.e. quantities that become conserved when the 
dissipative terms are taken to be exactly zero. For dissipative 
systems, an external force that acts as a source for the ideal in- 
variant quantity is employed and statistically stationary states 
of the system in regard to the ideal invariant quantity can be 
achieved. 

For hydrodynamic (HD) turbulence, the energy represents 
such an ideal invariant quantity along with kinematic helic- 
ity. In the inertial range, defined as the interval of scales 
smaller than that on which external forces act on and larger 
than scales where dissipative effects dominate, the properties 
of energy transfer functions, which depend exclusively on the 
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nonlinear dynamics of turbulence, are expected to have a gen- 
eral, universal behavior For a plasma medium, the nonlinear 
self-coupling problem of large spatio-temporal scales can be 
investigated in the magnetohydrodynamic (MHD) limit. For 
MHD turbulence, three quadratic ideal invariants exist: the 
total energy, cross-helicity and magnetic -helicity. Although 
strictly speaking only the total energy is an ideal invariant for 
MHD, it is often interesting to know the behavior of both the 
kinetic and magnetic energy transfer channels. 

In MHD turbulence, the locality problem is more complex. 
Not only multiple energy fluxes exist compared to just one 
for HD turbulence, but different dynamical states are now 
possible, providing different behaviors for the fluxes across 
scales. Traditionally, with each scale of motion £, which can 
denote the size of an eddy, a wavenumber k ^ 1/^ is as- 
sociated and the interactions between fields filtered to cor- 
respond to different scales Ufe(x,i) and hk{x,t) are inves- 
tigated. Alternatively, the same ideas can be applied to the 
Elsasser representation! 1 1 of the fields for which the nonlin- 
ear terms can be interpreted physically as the scattering of 
contra-propagating Alfven waves. In this approach, a Fourier 
space decomposition of the velocity field becomes the nat- 
ural framework and the transfer function between scales are 
obtained by selective integration (filtering) over the Fourier 
modes denoted by the wavevectors k. A review on the works 
employing this approach for MHD turbulence was done re- 
cently by Mininni[2|. 

Another way of looking at the problem is to relate the local- 
ity properties of the energy fluxes to the scaling of the velocity 
and magnetic fields. It is common[3,,4J to consider the order 
of the fluctuation for a scale Svk to depend on the energy E{k) 
as 5vk = [kE{k)]^^^ . Assuming a Kolmogorov type scal- 
ing for the energy spectrum in the inertial range, Exik) = 
C/f e^/^/c^^/"^, we obtain the scaling for the velocity field to 
be 5vk ~ k~^/^. Using the Iroshnikov-Kraichnan form of 
the energy spectrum characteristic of weakly Alfvenic turbu- 
lence, EiK{k) = CiK{VAef'^k-^/^, we find 5vk ^ fc"^/^ 
where Va — B / y/Airp is the Alfven speed of the guide mag- 
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netic field. Since these two scaling laws correspond to the two 
known limits of MHD turbulence tSJ, we can write in general 
that Svk ^ fc^S where ^ G [1/4, 1/3]. 

Using the smoothness condition for the MHD fields, which 
can be seen as arising from the scaling index <r, Aluie and 
Eyink showedJS) that for a dyadic (octave) separation of 
scales, the energy fluxes in MHD are local and can be charac- 
terized in the infrared (IR) and ultraviolet (UV) ranges by the 
scaling exponent with a limit depending on of the two fields. 
In particular for = 1 /3 the IR and UV limits are found to be 
±2/3. Since these limits are related to the MHD scaling, they 
are quite robust as results. 

In a work|7|, the present authors reported a locality expo- 
nent value of 1/3 for the energy conversion flux appearing in 
the kinetic energy equation and which that is responsible for 
the conversion of kinetic energy into magnetic energy. This 
limit would imply a stronger nonlocal effect for MHD tur- 
bulence. Since in the past it was suggested that this value 
could be a result of the external force polluting the inertial- 
inductive range, additional simulations have been performed 
and the values of the locality exponents are readdressed in this 
paper 

II. MHD STATIONARY STATES 

The locality properties of the scale fluxes are investigated 
for stationary state solutions of the incompressible MHD 
equations: 

= -u - Vu + b - Vb + i/V^u + f"- Vp, (1) 

ot 

— = -u • Vb + b • Vu + TyV^b + f , (2) 
ot 

V-u = 0,V-b = 0, (3) 

where u = u(x, t) is the fluid velocity field, b — b(x, t) is 
the magnetic field expressed in Alfven units and p — p(x, t) 
is the total (hydrodynamic + magnetic) pressure field divided 
by the constant mass density, p. Due to the incompressibility 
condition, the pressure p is an auxiliary variable and can be 
formally eliminated by solving the Poisson equation, 

V^p = -Vu : Vu + Vb : Vb . (4) 

Throughout this work, the fluid viscosity and the magnetic 
diffusivity 77 are taken to be equal. 

The zero divergent, external force fields f " = f "(x, t) and 
£6 = f''(x, t) act on the velocity and magnetic fields, respec- 
tively. The two forces are part of a forcing mechanism that im- 
poses the injection rates of the MHD ideal invariant quantities. 
We will refer to this type of forcing mechanism as a kinetic- 
inductive force. A kinetic only forcing method (f" = f and 
f = 0), used previously in the literature for similar studies 
ll8][9l, is also employed. In Fourier space, the forces f"(k) 
and {^(k) are assumed to be local, zero divergent quantities 
that act equally on all the modes within the wavenumber shell 
Sf = [kinf, fcsup]- Usually, the shell Sf is considered at large 
scales and is sufficiently thick to contain a large number of 



modes so that no preferential direction is introduced in the 
flow. In essence, both the kinetic only force f and the kinetic- 
inductive forcing mechanism f { " are considered to be pro- 
portional to the fields as, 

f (k) = Ci(k)u(k) + C2(k)cl>(k) , (5) 

where a; = V x u is the vorticity and 

f (k) = Cf"'^^ (k)u(k) + Cj"'''^ (k)b(k) . (6) 

The real valued parameters Ci and C2 are obtained by re- 
quiring the force to inject into the system a certain amount of 
energy in unit time (power) e and zero kinetic helicity. Simi- 
larly, for the kinetic-inductive version of the force the 

and cj" ''^ parameters are obtained by selecting the amount 
of energy and cross-helicity (ecr) injected into each of the two 
fields. The cross-helicity parameter a is bounded in the in- 
terval [—1,1]. Since the forces parameters are real, the forc- 
ing methods do not influence phases of the fields, which en- 
sures that no change is made in the type of turbulent structures 
present in the system. This might generate different results 
compared to the case of injecting cross-helicity by imposing 
the alignment of u and b in the real space for example. The 
advantage of the type of forces used in our work is in the abil- 
ity to select the levels at which the ideal invariant fluxes will 
relax to once the stationary state is achieved. Since the forces 
are proportional to the fields, the characteristic time of the 
force will tend to be equal to that of the nonlinear cascade for 
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FIG. 1. The compensated energy spectra for run I (a), run III (b) and 
run IV (c). 
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Run 


Nr. of modes 


^max 


V — r] 


Force Type 




e 


O" 


Rx 


rA 








I 


512^ 


256 


0.00050 




[2.5, 3.5] 


0.1 


0.0 


217 


0.47 


1.55 


1.86 


1.05 


II 


512^ 


256 


0.00050 


r&f'' 


[2.5, 3.5] 


0.1 


0.4 


296 


0.58 


1.92 


1.94 


1.06 


III 


512^ 


256 


0.00055 


f " only 


[1.5,3.1] 


0.32 


0.0 


362 


1.96 


1.76 


0.37 


1.22 


IV 


1024^ 


512 


0.00030 


f " only 


[1.5,3.1] 


0.32 


0.0 


450 


1.72 


1.94 


0.35 


1.15 



TABLE I. The simulations control parameters and the relevant diagnostics at the moment the flux analysis is being performed. The diagnostics 
are defined as: Rx = 2..™^ s/^J^^^y = / E\ L^^"^ = 2.1^^|^^^ and , = [.3/,] V4. 



a scale located in s/, independent of the energy injection level 
selected. A more detailed description of the electromagnetic 
forcing mechanism is given in another work[ 10|. 

The equations ([T]|2|, with the appropriate choice of forc- 
ing, are solved using a pseudo-spectral solver, using periodic 
boundaries conditions in all three directions. The size of the 
domain in each direction is 27r. The solver uses a FFT algo- 
rithm for the space discretisation and a order 3rd Williamson- 
Runge-Kuta method for the time advancement. The time step 
is computed automatically to be consistent with the CFL cri- 
terion. The nonlinear terms are partially dealiased using a 
phase-shift method [11]. The simulations are run until a sta- 
tistically stationary state is reached for both the velocity and 
the magnetic field. Multiple runs were made, details being 
listed in Table|l]for the well resolved stationary states investi- 
gated. Run IV is obtained from run III by increasing the nu- 
merical resolution and decreasing the viscosity and magnetic 
diffusively accordingly. The compensated spectra for the runs 
considered are displayed in Fig{T] Although all the data has 
been investigated, runs I and III will be mostly used for data 
display to exemplify best the effects generated by the two type 
of forces. 



field notations {X, Y, Z} stand in for m or & depending on the 
specific transfer and only their position in the transfer expres- 
sion is important. While Z represents the advecting fields and 
Y the advected field, X stands in for the receiving field. The 
delta Dirac function limits the transfers to wave-vector triads 
that satisfy the conditions k + p + q = 0. Note that in our 
previous work on this topic|7|, terms in the magnetic energy 
equation were denoted by Tut and T^u, i.e., with subscripts u 
and b interchanged compared with notation used in (j9]l. 
The equations (|7]|8]l also contain energy injection terms 
(k) which are due to the external forces and have the form: 

/^(k)-5R{f^(k).X(-k)} , (10) 

For stationary state turbulence, the total energy injection level 
imposes the total flux level in the inertial range and is given 
by the summation of the two contributions: /" + Individ- 
ually, /" and do not constrain respectively the the kinetic 
and magnetic energy flux levels, as the kinetic and magnetic 
energy are not themselves conserved quantities. For stationary 
state turbulence, the kinetic and magnetic energy fluxes relax 
to a level constrained by the kinetic and magnetic dissipation 
rates, respectively. 



III. LOCALITY FRAMEWORK 

A. Energy equation for a mode 

Since we are interested in the study of locality between 
scales, the IVIHD equations are solved in Fourier space. The 
energy equations for a mode k are easily derived and read as: 



dt 
dt 



- 7^:,u(k) - n[,{k) - 2,yk^E^k) + /"(k), (7) 
TL(^) - tU^) - 2^eE''ik) + , (8) 



Each of the non-linear terms of the type — Z • V Y appearing 
in the right hand side of the field X evolution equations, will 
generate energy transfers in spectral space of the form: 

?fz(k) = JJdqdp^\lk- Z(q)Y(p) ■ X(k)} X 

(5(k + p + q), (9) 

where Jft stands for the real part of a complex number and we 
have used the reality condition X*(k) = X(— k) which re- 
sults in Tyzi^) — Tyzi—^) for the real valued transfer. The 



B. Triad transfers and conservation properties 

From Eq. (j9]l, we see that the net energy received by the 
mode k of field X is due to the interaction with all possible 
modes p and q which form a triad. It is useful to look at the 
energy transfer for a individual triad, defined as: 

r^^(k|p, q) = ^^{zk ■ [Z(q)Y(p) + Z(p)Y(q)] • X(k)} , 

(11) 

for k + p + q = and zero otherwise. Since the triad trans- 
fer function is symmetric in p and q, i.e. TY^{k\p,q) = 
r^2(k|q, p), we have expressed this in an explicit way in 
the definition (111. Because of this symmetry, although we 



know the role played by each field in the interaction we have 
no way of differentiating between the contribution of modes p 
and q in a unique way. Following the work of Dar and Verma 
llT2HT3ll . a mode-to-mode transfer function can be introduced 
by accounting for a circulating transfer (an uncertainty) that 
cancel itself exactly for the triad transfer. 



5#,^(k|p|q) = 5R{[zk.Z(q)][Y(p) . X(k)]} 



(12) 
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for k + p + q = and zero otherwise. The triad transfer can 
be expressed now as the sum of two mode-to-mode transfers: 



Y,Z 



(k|p,q) 



1 



{^^^(k|p|q) 



^^z(k|q|p)} 



(13) 



For the energy transfer, although the nonlinear terms have the 
same form, they have different physical significance. The 
terms where u is the advective quantity, are respectively re- 
sponsible for the kinetic and magnetic energy conservation. 
This fact is expressed using the triad transfers as a conserva- 
tion of interaction between the modes that make up a triad: 



Tr,„(k|p,q) +T^,Jp|q,k) +T„"_„(q|p,k) = , (14) 



and 



Tfc'',„(k|p,q) + r,'',„(p|q,k) +r,^^Jq|p,k) = . (15) 

The other two terms, where h is the advective quantity, do not 
conserve the energy individually as they represent transfers 
from one field to the other (cross-field transfers). However, 
their sum, which accounts for the cross-field transfer, is con- 
served for any triad, 

T^bCklp, q) + r,';,(p|q, k) + r,':,(q|p, k) + 

T„%(k|p, q) + r„%(p|q, k) + r„^(q|p, k) = . (16) 

The cross-field transfers are responsible for the conversion of 
kinetic energy into magnetic one and vice-versa and it is due 
to their existence that the kinetic and magnetic energies are 
not conserved individually. This aspect needs to be considered 
when looking at the fluxes generated by the two terms taken 
separately. 

As a side note, we see that the same mechanism exists for 
cross-helicity, which has units of energy and up to a point can 
be considered as the energy contained by the velocity mag- 
netic interaction. For the cross-helicity, the transfers T^f,, 



T" J account for the transfer of information for the same field, 
while „ + T^^ represents the transfer between the veloc- 
ity and magnetic field. These cross-helicity interactions are 
conserved for any triads. For a triad, the four nonlinear terms 
appearing in the MHD equations (l]|2l generate eight trans- 
fers of the type (111, four in the energy equation and four in 
the cross-helicity equation, which contribute to six conserved 
interations, three in the energy equation and three in the cross- 
helicity equation. When using the Elsasser variables, the re- 
sulting transfer terms are just the combination of the eight 
transfers appearing for u and 6 representation and nothing 
more. If we start in Elsasser formalism and desire to recover 
fully the u and h representation, the residual energy (corre- 
lation of the co-propagating and contra-propagating Alfven 
waves phase velocities) needs to be taken into account. 



C. The scale transfer functions 

To quantify the transfer between scales, we decompose 
the wavenumber space into a series of disjoint shells Sfz = 



(/cx;_i, k]c\, similar to other works on the subject lfT4HT8l . 
The velocity field u'^ and magnetic field b'^ (note that the hat 
denotes Fourier transform) contained in a shell identified by 
the index k are found, by the use of a sharp spectral filter, to 
be: 




u(k) if |k| e stc 

if |k| i 

b(k) if |k| e s,c 

if |k| ^ 



(17) 
(18) 



The choice of a sharp filter compared to a smooth one has been 
investigated by Domaradzki and Carati lfT9l l20ll who found 
that the transfer functions and energy fluxes are similar in the 
two cases for sufficiently compact smooth filters. The shell 
wavenumber boundaries are obtained from the geometrical 
progression: k/c — 4: x 2^'^^^)/'*, with fco — 0. We obtain 
iV = 25 shells (/c = 1, 2, N) for /c„iax 256 and = 29 
for fcmax = 512. The use of a geometrical progression for 
the shells boundaries, assures us of capturing a sufficiently 
localized signal in both spectral and real space. For a unit 
linear separation for the shell boundaries, the angle integrated 
quantities would be recovered. When working with angle in- 
tegrated quantities, the designation band is usually employed 
instead of shells, as the wavenumber space decomposition is 
seen as selecting bands of the one-dimensional wavenumber 
space. The real space representation u'^(x) of the shell fil- 
tered velocity field, corresponds to a characteristic velocity 
for a scale Skjc, with the shell based scales Skjc separated as 
^kjc-i/Skic ^ 2^1^. This separation should be considered 
as the smallest separation between scales that we account for 
and not as a scaling directly linked to turbulence itself. In real 
space, the total field can be recovered by summing over each 
shell filtered contribution: 



N 



N 



u(x) = ^u'=(x), b(x) = ^b'=(x) 



(19) 



/c=i 



/C=l 



Numerically the transfer occurring between the shell fil- 
tered fields X'^, and Z^, is computed as. 



S^Mv\q) 



E 



sR{ik-Z2Y^(k)-X'=(k)}. (20) 



7C 



The triple transfer between shells 5y ^(/cIpIq), ignoring the 
fields nature but taking into account their position in the inter- 
action, has a more precise interpretation as the transfer to shell 
K from shell v, mediated by (through advection by) shell q. 
The contribution to a scale from the other two scales, regard- 
less on the role they play in the interaction can also be defined 
as. 



(21) 



and can be useful when employing infinitesimally thin shells, 
as the recovered symmetry would allow one to work with this 



quantity in direct analogy to the triad transfer (111. 
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The flux trough a sheU boundary (kc) is then defined by 
partial summing the transfer band spectra Tyzi^-)^ 

N N N N 

!C=c+l IC=c+lQ=lV = l 

(24) 

These fluxes for runs I and III are shown in Fig|2] 



D. The flux locality functions 

The infrared locality function is defined by taking a probe 
wavenumber boundary kp, so that kp < kc, and it measures 
the contribution to the flux trough kc from triads of modes 
with at least one wavenumber less than kp, 



N 



^iTY.zi^pl^c) ~ E/ 
/C=c+1 



■ N p 

E E S^,zMv\Q) 



N 



E E s^Mv\q) 

V=l Q=p+1 



(25) 



FIG. 2. The total flux for run I (a), run III (b) and run IV (c) and the 
contributions made by the four nonlinear tenns, displayed in absolute 
values. The vertical lines depict the shell boundaries. For run I, up 
to shell 7 and for run IV, between shell 6 and shell 12, the flux 11"^ 
is negative. 



The function Sy zM^^Iq) is computed numerically from 
solutions of the MHD equations and forms the basis of our 
analysis. Knowing Sy zM^^Iq) allows us to compute the 
shell-to-shell transfers by summing over all possible advec- 
tive shells. 



E^ 

c 



s) 



(22) 



The shell-to-shell function Py2;(ic\'p) can be seen as a trans- 
fer from shell v of field Y to shell k of field X and possesses 
the antisymmetry property Pyzitclv) = — P^2(p|/c). 

Similarly, the net transfer to a shell Tyz (ic) can be found 
by summation over v and q. 



w - E w^) = E E 

■P V Q 



(23) 



and can be seen as the shell integrated transfer spectra. When 
summing the net transfer ( [23] l over k, which is equivalent to 
integrating the nonlinear transfer over the entire space, we ob- 
tain zero only for the interactions that are conserved in a triad. 
This fact requires extra care from us when defining and inter- 
preting the energy fluxes through a shell boundary surface, as 
not all fluxes go to zero in the UV limit (large wavenumber 
limit). 



In the second term, the sum over shell q starts from p + 1 to 
avoid double counting. For the limit kp -> kc, we recover 
the flux for all terms. It is customary to normalize the local- 
ity functions to the flux trough kc, which from the definition 
shows that is one for kp — kc and decreases for kp/kc < 1, 



lim 



Hy z{kc) 



lim 



Ilii- Y zi^clkc) 



1 . 

(26) 



Similarly, the ultra-violet locality function nuv(fcp|fcc), which 
accounts for the contributions of small scales to the transfer 
can be defined. The UV functions are strongly influenced 
by the dissipation (or more correctly by the shape of the flux 
which naturally decreases for high k in dissipative systems) 
and do not exhibit a clear scaling. From our numeric analysis, 
a lower Umit of 1 /3 can be safely inferred for the UV local- 
ity exponent, without excluding the possibility that the actual 
asymptotic values are closer to 2/3. The same dependency of 
the scaling slope clarity on the shape of the fluxes is experi- 
enced by the IR locality function, as well. However, due to the 
forcing mechanism used, these effects are found to influence 
the IR locality exponent scaling to a lesser extent and for this 
reason we focus in this work on the analysis of the IR locality 
exponents. 

For the kinetic forced run, the IR locality functions are 
shown for the four energy fluxes in Figj3] For the flux 
of kinetic energy (11^ and the flux of magnetic energy 
(n^ jj) a 2/3 scaling can be observed. The conversion terms, 
i.e. kinetic to magnetic (HJ^^) and magnetic to kinetic (H'^^ j) 
fluxes, exhibit a more complex picture. For these two terms 
the locality scaling seems to be between 1/3 and 4/3. Look- 
ing now at the kinetic-inductive forced run, Fig|4j we see that. 
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while the kinetic and magnetic energy fluxes asymptote to the 
same 2/3 value, the conversion terms show a clear 1/3 scal- 
ing. Since only absolute values are plotted in these figures, the 
presence of a cusp is indicative of a change in the sign of the 
locality functions. As the locality functions are normalized by 
their respective flux through fee, the change in sign indicates a 
change in the nature of the contributions to the flux across the 
cutoff scale. 

To obtain a better physical understanding of the 1/3 scal- 
ing value, we plot in Fig|5]the IR locality scaling for the entire 
conversion term (containing the interactions of both of the two 
cross-field terms). We remind the reader that only when con- 
sidered in its totality does the cross-field interaction is energy 
conserving. Because of this conservation and cancelation ef- 
fects, the global scaling for the conversion flux drops to values 
closer to 2/3 or even larger denoting a more local behavior. 

IV. A DETAILED ANALYSIS OF AT THE TOTAL FLUX 

The cancelation nature of the fluxes and the effect this has 
on the locality properties require an additional study. To max- 
imize the lesson learned from HD turbulence and use the in- 



tuition gained by the community in this field, we look at the 
total flux. For MHD, we define the total net energy transfer to 
a shell as: 

T{K) = Tl^(K) + Tl^{K) ~ TI,{k) - T^K) . (27) 

For MHD turbulence, the total net transfer T{k) and the af- 
ferent quantities (e.g. total triple transfer between shells, total 
shell-to-shell transfer) are the ones that have a behavior sim- 
ilar to the respective quantities found for HD turbulence. For 
each of the quantities introduced in the previous section, a 
similar sum can be performed. Alternatively, one could just 
define the total energy triad transfer as the sum of the four en- 
ergy triad transfers possible and build up the framework for 
the total energy exchange among scales. From this point on, 
quantities for which we omit the X, Y, Z labels refer to total 
quantities. 

To better understand the contribution to a flux trough a shell 
boundary (fee), we look at the total flux, defined as: 

N N N N 
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For the three runs, the total fluxes are shown in Fig.|2] Writing 
in detail the second two sums in regard to the shell index c we 
find four terms, each one with a different physical interpreta- 
tion. 



N N 



V=l Q=l V = l Q = l 

c N 

+ E E }(") 

•P=l Q=c+1 
N c 

•P=c+1 Q=l 

N N 

+ E E ^(^l^ls) ' (29) 

•P=c+1 Q=c+1 

The first term [i) contributes to the flux as the transfer from 
large scales to small scales, advected by the large scales, while 
the second term [ii) represents the transfer from large scales 
to small scales, advected by the small scales. The third term 
[Hi) can be seen as the transfer caused by the large scales 
advection of the small scales and in general is found to be 



close to zero. The last term {iv), which could be seen as the 
transfer from small scales to small scales by the advection of 
the small scales, is exactly zero due to the conservation of 
the total energy for any closed set of modes interacting only 
among themselves. The contributions to the total flux can be 
seen in Fig|6] 

To account for the degree of locality of the total flux, we 
look at the locality properties of the two terms (fluxes) that 
bring the main contribution [i) and (ii), 

n(fc,)«n(')(fc,) + n(")(fce) , 

N c c 

n^'Hfcc)- E EE^W^Is)' 

/C=c+1 V=l Q=l 
N c N 

n^''Hfcc)= E E E ^(^l^ls)- (30) 

K=c+l V=l Q=c+1 

Knowing the value of the flux through a shell surface kc we 
want to know how much of this flux is due to modes with 
wavenumbers close but smaller then kc- For this purpose we 
take a probe (test) wavenumber boundary kp, so that kp < kc, 
and we measure the contribution to the flux trough kc from 
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FIG. 5. Infrared locality functions for the total conversion term 
Hir + Hii- t,i,(fc|fec); run I (a) and run III (b). 



modes with wavenumber less than fcp, similar to the philos- 
ophy of the IR locality functions. By keeping fee fixed and 
varying fcp, we should obtain a smaller and smaller contribu- 
tions. The rate at which the contributions become smaller is 
related to the locality of the flux. 

In the case of 11'^*' we can take the probe on the giver shell 
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FIG. 6. Total flux for run III (gray) and run IV, displaying the 
contributions made by the main terms. The iiii) terms are 10^'' 
times smaller compared to the other two nonzero terms and fluctuate 
around zero. 



FIG. 7. Locality scaling for the contribution of energy across the 
fcc, advected by the entire large scale region (Hp' (fcpjfcc)): run I (a) 
and run III (b). Inset pictures show the same pictures without the 
contribution of the first shell responsible for the forcing. 



index v, 

N p c 

n^;\k,\k,)^ E EE^w^is)' (31) 

K=c+lV=l Q=l 

or on the advecting shell index q 

N c. p 

n^^(fcp|fce)= E EE^W^Is)- (32) 

K = c+lV = l Q = l 

The locaUty nature of the partial flux II'^*-' is related to both 
effects, however, performing the probe variation separately al- 
lows us to understand the locality nature of these two effects 
individually. First, by varying the probe on the giver shell 
and keeping the advection to the entire large scales region, 
we transfer energy from increasingly separated scales. From 
Fig|7] we see that the locality function scales as 7/3, which 
would indicate a highly local behavior. Second, we vary the 
probe on the advecting shell while keeping the entire large 
scale interval as the giver of energy. In Fig|8]we see that the 
advection brings a more nonlocal behavior to the partial flux. 
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When considered together it is the more nonlocal behavior 
that dominates the asymptotic behavior. 

For the partial flux H^*'), a more pronounce nonlocal be- 
havior has been found. The scaling is close to a 2/3 value for 
the kinetic forcing and a 3/4 value for the kinetic-inductive 
forcing. The value of 3/4 for the IR locality exponent is con- 
sistent with the theoretical value found by Aluie and EyinkfS) 
for fields that scale as k~^l^. 

From the locality picture of the partial fluxes we see that 
the advecting effects are the ones responsible with the value 
of the IR scaling exponent for the total flux. Combining all 
these behaviors (advecting and advected) results in a scaling 
close to 2/3 obtained for Kraichnan IR locality functions for 
the total energy flux of MHD turbulence. 



V. CONCLUSIONS AND DISCUSSION 

After investigating the locality properties of MHD energy 
fluxes, it is found that asymptotically the dynamics tend to 
be dominated by local interactions. The nonlocal interactions 
that apparently exist, cancel themselves out. However, the 
locality is much weaker compared to the case of HD turbu- 
lence, which is characterized by the scaling exponent of 4/3. 
When using a velocity proportional force, two distinct expo- 
nents are observed for MHD turbulence, 1/3 and 2/3 for var- 
ious fluxes. The 1/3 exponent is even more obvious for the 
kinetic-inductive forcing mechanism. 

To better understand these values we investigated the to- 
tal energy flux. The lesson taught by the analysis of the total 
flux showed us that the locality of a flux can be seen as the 
locality of the contribution of each effect that makes up that 
flux. For similar level contributions, it is the most non-local 
channel in a flux that imprints the overall locality behavior 
For MHD turbulence, the overall locality seams to be close 
to 2/3 and confirms the analysis done by Aluie and Eyink[6^ 
who showed a 2/3 scaling behavior for the individual interac- 
tions of HD and MHD turbulence. Since for MHD turbulence. 



the interplay between velocity and magnetic field gives rise to 
much smoother correlated fields, the decorrelation effects due 
to the averaging procedure are less pronounced, which results 
in a smaller global locality exponent and justify the 2/3 value 
found for the various energy fluxes compared to the 4/3 value 
found for the energy flux of HD turbulence. This can be seen 
best by looking at the total flux, the corresponding flux for the 
velocity flux in hydrodynamical case. Various non-locality 
contributions and force influence cancel themselves, making 
the total flux scaling exponent the most reliable observation. 

The observed 1/3 value was initially thought to be caused 
by the forcing. Assuming that the 1/3 locality exponent is in- 
deed due to the external force or due to the pollution of the 
inertial range by the large scale strains and that for a proper 
inertial-inductive range the 2 /3 exponent index is recovered, 
this still represents a major problem for adequate modeling. 
A proper inertial-inductive range is so slow to appear that a 
huge number of modes need to be solved for any LES type 
model. For practical reasons, if the forcing range pollutes the 
inertial-inductive range to this extent it is better to take into 
account the more pronounced nonlocal behavior of the indi- 
vidual fluxes rather then trying to reach the inertial-inductive 
range. We believe that the 1 /3 scaling exponents is found due 
to the consideration of the conversion terms individually and 
is not just an artifact of the force. When triads are summed 
in a non-conservative way, fluxes of energy naturally appear 
This fact may represent another problem for practical mod- 
elling. In MHD shell modeIs ll2ni22l the two conversion terms 
are still accounted for separately (a stronger nonlocal charac- 
ter) but their interaction is still considered to be mainly with 
only the closest neighbor shells, in the same way as in the 
case of the more local kinetic and magnetic terms. Allowing 
the cross-field terms to couple to more distant shells might be 
more appropriate as the level of contribution neglected to a 
flux would be the same for all individual terms. 

A question remains: Is there any physical significance 
for the 1/3 exponent seen for the energy conversion terms? 
If we try to interpret MHD turbulence as scattering of 
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contra-propagating Alfven waves, then we actually work with 
the wave phase velocities regardless of actually employing 
Elsasser formalism or not (which just acknowledges the fact 
that we are mixing u and h in the definition of the wave phase 
velocities). For this interpretation there is no 1/3 scaling. 
Moreover, in this interpretation all the four energy fluxes and 
any combination of them are just partial fluxes, pieces that 
make up the entire physical picture. Only the energy of a wave 
should be considered in this case and that energy is the El- 
sasser pseudo-energy (which is equal to the total energy only 
in the absence of cross-helicity) and for which the respective 
fluxes have a scaling exponent close to 2/3. Complementary, 
if we just look at the interactions in a triad and assign a phys- 
ical interpretation to each conserved interaction {e.g. kinetic 
energy exchange, magnetic energy exchange, energy conver- 
sion), we need to consider the two conversion terms together 
and we again do not obtain the 1/3 scaling since the nonlocal 
contribution cancel themselves. 

However, when MHD equations for u and h are used, as is 
the case most often, the scaling for individual terms appears to 
be a legitimate question. Looking at the two energy evolution 
equations separately (notwithstanding the arguments that the 



two are not strictly redistributive by themselves) we can state 
that the kinetic energy of the system is modified by a flux 
(generated by the the entire nonlinear rhs) which possesses a 
1 /3 locality scaling exponent, since the u to 6 flux is larger in 
magnitude and has the strongest nonlocal behavior This in- 
formation may not be adequate for the physical interpretation 
of the energy transfers, but it may be useful for modeling if u 
and h equations are used and one tries to model effects of each 
term separately. 
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